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We investigate electronic and optical properties of the V-based ladder compounds NaV2 0s, the 
iso-structural CaV20s as well as MgV20s, which differs from NaV20s and CaV^Os in the c axis 
stacking. We calculate ab initio the A 9 phonon modes in these compounds as a basis for the 
investigation of the electron-phonon and spin-phonon coupling. The phonon modes together with 
the dielectric tensors as a function of the corresponding ion displacements are the starting point for 
the calculation of the A 9 Raman scattering. 

PACS numbers: 71.15.Mb, 71.27.+a, 63.20.K 
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I. INTRODUCTION 

The vanadium ladder structures NaV^Os, CaV^Os, 
and MgV205 are low-dimensional compounds with un- 
usual physical properties due to a strong coupling of 
charge, spin, and lattice degrees of freedom. NaV205 and 
CaV205 are iso-structural (space group (Pmmn)) 
while MgV205 has the symmetry group Cmcm. The main 
structural element of all these compounds is a simple lad- 
der containing V-0 rungs and V-0 legs providing a possi- 
bility for electrons to hop within a rung and between the 
rungs of one ladder, whereas the hopping matrix element 
between different ladders is considerably smaller. The 
highest occupied states are formed by d xy electrons of V 
mixed with oxygen 2p orbitals. CaV^Os and MgV205 are 
both half-filled compounds with similar electronic prop- 
erties, i.e. one d xy electron per V site. However, they 
show different magnetic behavioni. The quarter-filled 
NaV205 where one d xy is shared by two V sites within a 
rung 2 , undergoes a spin-Pcicrls transition at T c — 35 K 3 
accompanied by a charge ordering and a large lattice dis- 
tortion. In this compound, phonons induce strong charge 
fluctuations near T c that lead to fluctuations in the spin- 
spin exchange J^. Therefore, investigation of lattice dy- 
namics and electron (spin) - phonon coupling in ladder 
compounds can provide a clue to their properties. 

Investigations of these compounds by optical spectro- 
scopies - ellipsometry and Raman scattering - have shown 
many interesting features. The former gives the dielectric 
function, whereas the latter provides information about 
phonons and low-energy excitations and their coupling to 
light. For NaV2C>5 the phonon Raman spectra demon- 
strate very large changes below the transition points such 
as variations in the electronic background^ and the ap- 
pearance of intense new Raman active modes indicat- 
ing a strong coupling of the lattice to the order param- 
eter of the transition 4,5 . With the aim to understand 
similarities and differences in the optical properties of 
NaV205, CaV205, and MgV^Os we perform ab initio 
investigations within density functional theory (DFT). 
We calculate the fully symmetric A g phonon modes and 



their influence on the electronic structure and the di- 
electric function. Such phonon-induced modulations of 
the polarizability are a measure for electron-phonon cou- 
pling and the source of phonon Raman scattering of light. 
Although the vanadates are considered to be correlated 
materials^, the results obtained in the framework of DFT 
show good agreement with experiment in terms of their 
lattice dynamics and optical properties. Since similar 
findings have been reported earlier for high temperature 
superconductors 8 such investigations are justified. 



II. COMPUTATIONAL DETAILS 

All calculations are performed within DFT using the 
full-potential augmented plane waves + local orbitals 
(FP-APW+lo)S formalism implemented in the WIEN2k 
codeiS. Exchange and correlation terms are described 
within the generalized gradient approximation (GGA)^ 1 ". 
All atomic positions of the three structures have been 
relaxed as a starting point for the investigation of the 
phonon modes within the frozen-phonon approach. We 
have displaced the atomic coordinates from their equilib- 
rium positions, where for each degree of freedom two dis- 
placements in both, positive and negative direction, re- 
spectively, have been used. The corresponding forces es- 
tablished the dynamical matrix for the A g phonons. The 
dielectric tensors are computed within the random phase 
approximation with the Kohn-Sham eigenstates taken as 
an approximation to the single particle states. 



III. DIELECTRIC FUNCTION 



Fig. n and Fig. [21 demonstrate the diagonal in-plane 
components of the dielectric function £ xx (lu) and s yy (uj) 
of the three compounds investigated. The real part of 
s(u)) was obtained by the Kramers - Kronig relation from 
its imaginary part, where no broadening of the optical 
transitions was assumed. The strong absorption peak in 
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FIG. 1: The dielectric function e^fw) for NaV 2 5 , CaV 2 5 , 
and MgV205. The :r-axis is parallel to the rungs. 



FIG. 2: The dielectric function £ yy (uj) for NaV 2 5 , CaV 2 5 , 
and MgV 2 05. The j/-axis is parallel to the ladders. 



the imaginary part of e xx (uj) at u> close to 1 eV and an- 
other weak feature in Im e yy {ui) at u> ~ 1.25 eV are due to 
transitions between the bands originating from bonding 
and antibonding states within one rung separated by « 
2|tj_| in energy, where t± is the in-rung hopping matrix 
element. Due to the dispersion of the bands the peak 
position is larger than 2\t±\. Our results correspond to 
t± = —0.35 eV, in agreement with Ref. |2|. At the same 
time, a peak at 1.75 eV in Im s yv {ui) arises from tran- 
sitions within the legs. For NaV20s the calculated di- 
electric functions are in good agreement with available 
experimental data™. CaV205 demonstrates similar fea- 
tures, while for MgV205 the interpretation is more com- 
plicated. 



IV. A 9 PHONON MODES 

The A g phonon modes preserve the symmetry of the 
unit cell and in particular the inversion symmetry, and, 
therefore, are seen in the Raman spectra. Table |U 
presents the calculated and measured eigenfrequencies of 
these modes for all three compounds. In all modes both, 
V and O displacements are considerably involved in the 
eigenvectors. The apical oxygen ions only contribute to 
the highest-frequency modes (w p h > 900 cm -1 ), thereby 
vibrating along the z-axis. The theoretical eigenfrequen- 



cies of NaV2C>5 and CaV20s for all but the lowest mode 
are in very good agreement with experiment, while the 
agreement for the MgV20s modes is worse. The eigen- 
vectors of all three compounds correspond well to the 
assignment given in Ref. Il3l 



PHONON-INDUCED CHANGES IN THE 
DIELECTRIC FUNCTION 




FIG. 3: Illustration of the in-plane 02-V stretching phonon 
mode. The lengths of the arrows are proportional to the re- 
spective components of the eigenvector. 
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TABLE I: Experimental and calculated phonon frequencies uj p h for all three ladder compounds and the assignment of the 
modes. 
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FIG. 4: Effect of the V-02 stretching on the dielectric func- 
tion in NaV 2 Os. 



FIG. 5: Effect of the V-02 stretching on the dielectric func- 
tion in CaV 2 Os. 



Raman scattering arises due to phonon-induced 
changes in the band structure, and, correspondingly, in 
the dielectric function 5e aa (uj), where a are the Carte- 
sian indices. The intensity of the one-phonon scattering 



Ir for the incident and scattered light parallel to the a 
axis can be written as 



Ir 



:n^<5io> 



in 



4 



Here Q is the coordinate of the phonon mode, and the 
matrix element represents transitions between the states 
|0) and (1|, which are the phononless ground state and 
the one-phonon state, respectively. 

Fig. |31 represents the ion displacements in the V-02 
stretching mode which is taken as an example to demon- 
strate the influence of the phonons on the dielectric func- 
tion. Figs. 01 and show the changes in the xx and 
yy components of the dielectric function of NaV20s and 
CaV2C>5, respectively. The black lines correspond to the 
situation of the undisplaced positions, whereas the grey 
lines reflect the distorted lattice with the displacement of 
the 02 ion along the x-axis taken to be 0.02 A, which is of 
the order of magnitude of the zero-point vibrations. The 
V ion is shifted correspondingly according to the eigen- 
vector of the mode. As it is seen in the figure, this mode 
causes big changes in s xx (u) at u> close to the resonance 
at 1 eV which corresponds to the intra-rung transitions 
due to a modulation of the t± hopping within a rung. 
A general feature of the phonon-induced changes in the 
dielectric functions is that all the peaks are considerably 
shifted. This fact implies large values for the first deriva- 
tives of the dielectric functions, which govern the Raman 
scattering intensities. We therefore predict strong reso- 
nance effects in phonon Raman scattering close to these 
peak positions. E.g. the 1 eV region of NaV20s should 
exhibit resonances in the xx component. For CaV20s we 
expect even more resonances since the peak structure of 
the dielectric function is more pronounced. Due to our 



findings we suggest Raman scattering experiments to be 
performed below 2 eV, the energy region where we expect 
very high Raman intensities. At the same time, we can 
explain the experimentally observed intensity enhance- 
ment for NaV2C>5 when changing the frequency of the 
laser beam from 1.95 to 2.4 eV in xx polarization^, since 
the latter frequency is closer to the peak in Re e xx (co) at 
~3 eV as can be seen in the upper panel of Fig^ 



VI. CONCLUSIONS 

We have performed ab initio calculations for the opti- 
cal properties and A g phonon modes in the ladder com- 
pounds NaV 2 5 , CaV 2 5 , and MgV 2 5 . For NaV 2 5 
the calculated e(u) compares well with available experi- 
mental data. For NaV2C>5 and CaV20s the phonon fre- 
quencies are in a very good agreement with the experi- 
ment, while for MgV205 the agreement is not as good. 
Due to the strong influence of the ion displacement on 
the dielectric function we predict pronounced resonance 
effects in phonon Raman scattering in the energy range 
below 2 eV. 
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